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Comparison  of  theoretical  single-wire  current  response  over  an 
ideal  ground  plane  and  response  of  a single  wire  over  treated 
sand  ground  plane  


This  is  the  third  in  a series  of  reports  treating  the  multiwire 
array  and  designed  to  meet  the  needs  of  an  EMP  vulnerability  assessment 
and  hardening  program.1'2  The  approach  to  the  shielded  multiwire  cable 
is  extended  to  the  wire  array  over  a ground  plane.  A limited  code  is 
presented  that  calculates  the  wire  current  distribution,  degree  of  wire 
coupling,  and  the  effective  single  conductor  equivalent  for  an  array 
exclusive  of  end  loading  effects.  Theory  is  compared  with  experiment 
for  the  case  of  wire  arrays  near  a less  than  ideal  ground  plane. 


2.  THEORY 


For  a parallel-wire  array  over  a lossless  ground  plane  (and  this 
treatment  disregards  wire  end  loading) , the  potential  at  any  wire  can  be 
written  in  the  form 


where  I^x.t)  is  the  current  on  wire  M and  M' 
With  the  ground  surface  forming  the  image  plane 


is  the  total  wire  count 
, Zjj.,  is  (in  mks  units)1 


1/2* cl/ ztnc 


1 Michael  J.  Vrdbel , The  Multiuire  Shielded  Array — Theory  and  Code , 
Harry  Diamond  Laboratories,  HDL-TR-1873  (October  1978). 

2Michael  J.  Vrabel,  Ueer'e  Manual  for  SCWAR — A Multiuire  Shielded 
Cable  Code  for  a System  Oriented  EMP  Vulnerability  Aeeeeement  Program, 
Harry  Diamond  Laboratories , HDL-TR-1887  (May  1979). 


radius  of  conductor  N) 


velocity  of  light  in  a vacuum 


This  treatment  is  restricted  to  the  case  of  a uniform 
constant  for  the  mediun  above  the  ground  plane. 


where 


K - a coupling  term  whose  exact  form 
treatment,  and 


this 


the  height  of  wire  N above  the  ground  plane 


Equations  (1),  (2),  and  (3)  define  a set  of  M'  simultaneous 
equations  that  can  be  solved  for  the  relative  current  on  all  array 
wires. 


The  impedance  to  ground  of  a single  wire  of  the  array  is  given  by 


The  effective  impedance  of  the  wire  array  to  ground,  is  given  by 


When  the  array  is  over  an  earth  ground  plane,  two  significant 
difficulties  arise: 

1.  The  concept  of  transmission-line  termination  for  a very  lossy 
system  is  ill-defined.  Since  the  development  of  a wire-array 
code  along  the  line  indicated  requires  a calculation  of  an 
effective  single-wire  response  with  no  reflection  off  the 
terminations,  error  is  introduced  with  this  calculation. 

2.  In  general,  the  losses  associated  with  an  earth  ground  plane 
play  a dominant  role  in  establishing  the  late-time  response  of 
the  array.  Since  a time  domain  code  can  handle  such  losses  only 
approximately,  the  array  late-time  response  will  be  distorted. 

Within  the  context  of  developing  an  BMP  systems-related  code,  the 
present  approach  can  be  justified  on  several  grounds. 

1.  It  permits  the  development  of  a reasonable  code  to  handle  a 
potentially  very  complex  subject. 

2.  The  uncertainties  introduced  by  the  late- time  response  are 
dominated  by  those  resulting  from  other  aspects  of  a 
vulnerability  assessment  program. 

3.  The  early-time  response,  the  region  of  most  interest,  is 
unaffected  by  loss  term  uncertainties^  introduced  with^  this 
approach  but  rather  reflects  the  accuracy  (or  lack  thereof^  of 
the  single-wire  code. 

This  treatment  assumes  that  the  cross-section  dimensions  of  the 
array  are  fixed  and  that  equilibrium  conditions  exist  for  all  wires. 
The  proper  measure  of  the  latter  assumption  is  not  the  rise  time  of  the 
incident  signal  but  rather  the  wire  responses  as  measured  at  either  end 
or,  because  all  coupling  relates  to  the  response  of  an  equivalent  single 
wire,  the  response  of  the  equivalent  conductor. 


3.  A COMPARISON  OF  THEORY  WITH  EXPERIMENT 

Experiments  were  designed  to  test  the  predictions  of  section  2.  The 
experiments  were  performed  at  the  Harry  Diamond  laboratories  scale  model 
facility  employing  a lossy  ground  plane  composed  of  sand  treated  with  a 
solution  of  sodium  and  calcium  chloride.  Soil  resistivity  measurements 
were  not  available.  The  effect  of  the  poor  ground  conductivity  is 
demonstrated  in  figure  1,  which  shows  the  experimental  and  theoretical 
response  of  a single  wire  as  a function  of  height  over  a ground  plane, 
all  data  normalized  to  the  maximum  height  of  measurement.  For  this  and 
all  subsequent  measurements,  the  signal  source  is  a horizontal  dipole 
near  the  ground  plane  (angle  of  illumination  10  to  20  deg)  at  a range  of 


vm* 

V * , * \» 


approximately  2 m.  The  measured  response  is  the  coupled  early-time  peak 
current  due  to  a broadside  illumination  with  a 5-ns  square  pulse*  The 
rise  time  of  the  illuminating  signal  was  such  that  equilibrium 
conditions  were  maintained  across  the  wire  array. 
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Figure  1.  Comparison  of  theoretical  single-wire  current  response  (solid 

curve)  over  an  ideal  ground  plane  and  response  of  a single  wire  over 

treated  sand  ground  plane  (both  curves  normalized  to  the  45.7-cm 

**.  • *•  -*  • • - .......  _»  „ ....... 

height) . 


Figure  2 shows  the  results  of  a series  of  two-wire  measurements. 
The  plotted  data  represent  the  ratio  of  the  peak  response  of  a single 
wire  to  the  response  of  the  same  wire  in  the  presence  of  a second 
wire.  The  abscissa  is  the  distance  between  the  two  wires.  For  all 
points,  the  measuring  wire  was  unmoved  and  the  second  wire  was  varied  in 
a vertical  plane  both  above  and  below  the  measuring  wire  and  in  a 
horizontal  plane  both  in  front  of  the  measuring  wire  (that  is,  between 
the  measuring  wire  and  the  dipole  antenna)  and  behind  it.  The 
measuring  wire  was  fixed  10.2  an  above  the  ground — a region,  as  shorn 
in  figure  1,  where  the  effect  of  finite  soil  conductivity  is  becoming 
significant. 


A series  of  three-wire  measurements  was  performed.  The  geometry 
employed  and  the  response  data  are  given  in  figure  3.  The  response  was 
measured  as  with  figure  2,  but  with  two  rather  than  one  additional 
wire.  The  open  circle  indicates  the  wire  upon  which  all  measurements 
were  performed. 


Figure  2.  Ratio  of  current  response  (S^  of  single  wire  in  two-wire 
array  to  current  response  (S2)  of  single  isolated  wire  (measuring  wire 
fixed  at  10.2  an  above  ground  plane).  Curve  A (theoretical)  and  points 
+ (experimental) — second  array  wire  varied  in  a vertical  plane  above 
measuring  wire.  Curve  B (theoretical)  and  points  o and  • 
(experimental) — second  array  wire  varied  in  a horizontal  plane  in  front 
(•)  of  measuring  wire  and  behind  (o) . Curve  C (theoretical)  and  points 
A (experimental) — second  array  wire  varied  in  a vertical  plane  below 
measuring  wire. 

*,/*»  »,/*,  wet  GEOMETRY 

EXPEMMEMTAL  THEORETICAL 


Figure  3.  Experimental  and  theoretical  current  response  of  a three-wire 
array.  S^:  response  of  a single  wire  in  the  array;  S2:  response  of 

same  wire  with  remaining  two  removed.  Drawings  are  scaled 
representation  of  wire  arrays  (wire  diameters  not  to  scale);  open 
circle  is  the  measuring  wire— fixed  for  all  measurements  at  10.2  an 
above  ground  plane. 
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4.  A CODE  FOR  AN  INFINITE  WIRE  ARRAY  OVER  A GROUND  PLANE 

A FORTRAN  code  was  developed  to  implement  various  provisions  of  this 
development.  Among  its  uses  are  the  following. 

1.  In  any  parallel  array  of  wires  it  must  be  ascertained  whether  a 
single-wire  or  a multiwire  array  code  is  required.  This  is 
determined  by  calculating  the  ratio 

*N  = ZNNIN A ZNMV  M = '•  2'  * * ‘ M*  <6> 

M 

for  all  wires.  For  Rjj  = 1,  there  is  no  coupling  between  wire  N 
and  the  rest  of  the  array.  The  smaller  this  ratio  becomes,  the 
greater  the  degree  of  coupling.  It  remains  to  the  user  to 
select  the  point  where  the  single-wire  treatment  is  no  longer 
valid  (this  judgment  being  based  on  overall  system 
requirements).  One  caveat:  the  calculation  is  based  on  a good 
ground  plane;  for  a lossy  system,  the  level  of  wire  interaction 
becomes  larger  as  the  wire  array  is  positioned  closer  to  the 
ground  plane  than  is  predicted  by  R^.  A more  stringent 
requirement  for  RN  is  required  for  an  array  near  the  earth  than 
for  an  aerial  array. 

2.  For  special  array  configurations,  where  the  wire  impedances  and 
end  loads  are  similar,  a circuit  analysis  can  be  performed  on 
the  array-equivalent  single-conductor  synthesis.  An  example  of 
such  an  analysis  appears  in  an  Army  Corps  of  Engineers  study.3 

3.  The  code  outputs  the  relative  current  distribution  among  the 
wires  for  the  infinite  array  case.  While  this  fact  has  some 
intrinsic  interest,  caution  must  be  exercised  in  interpreting 
these  data  since,  for  most  cases  of  interest,  wire-end  loading 
plays  a dominant  role  in  determining  wire  signal  levels. 

The  infinite  parallel-wire  array  code  outputs 

1.  all  input  data, 

2.  the  level  of  wire  interaction, 

3.  the  relative  current  distribution  among  the  array  wires  for  the 
infinite  transmission  line, 

3EMP  Disconnect  and  Protectiori  for  Cormercial  Power  Feeders— Otie  APB 
Pave  Paws  Radar,  U.S.  Army  Corps  of  Engineers,  HNDSP-  77-35 8 -ED- SR  (30 
September  1977). 
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4.  the  impedance  to  ground  for  all  array  wires  for  the  infinite 
transmission  line,  and 


5.  the  impedance  to  ground  of  the  wire  array  for  the  infinite 
transmission  line. 


The  code  input  data  are 


RAD(N) 


height  of  wire  N above  the  ground  plane,  and 


horizontal  location  of  wire  N along  the  ground  plane 
measured  with  respect  to  any  arbitrary  but  fixed 
reference  point.  X(N)  must  be  a positive  value. 


The  code  is  designed  to  handle  a maximum  of  30  wires  although  this 
limit  can  be  increased  by  altering  the  dimension  statements. 


Any  consistent  set  of  units  may  be  used  for  all  input  data.  A 
listing  of  the  FORTRAN  program  follows.  Near  statement  number  6 
subroutine  GELG  is  called;  this  is  an  IBM-supplied  subroutine  for  the 
solution  of  a set  of  simultaneous  linear  equations. 


DIMENSION  RAD  1 30 1 • 
DIMENSH  ! COR  130), 
DIMENSION  VOLT  C 30  I 


r(30l.  1130, 301,  COE Fl (9001 , C0EFR(30) 
ACURI30  I 


NAME  LI  ST /LI  ST A/NUN, A AD.X  ,Y 


C NOH  IS  THE  TOTAL  NUMBER  OF  CONDUCTORS 


X(M),  YIMI  ARE  THE  POSITION  VECTORS  OF  THE  HIRES  MEASURED  MITH 
RESPECT  TO  THE  GROUND  PLANE  AND  HIRE  1 


HR  IT  E 16 ,22  I 

FORMAT!*  MULTIPLE  PARALLEL  LINES  ABOVE  GROUND 
HR  IT  E (6 ,23  I 
HR ITE (6,18 INUM 

FORMAT! I0X  ,20HNUMBER  UF  HIRES  ARE  ,12/1 

00  26  N *1 ,NUM 

HRITEI6,19)N,RAD(N) 

FORMATM0X.1SHRADIUS  CF  HIRE  , 1 2 ,4H  is  .E10.4I 


VARIOUS  PARM.*| 


COM  INUE 
HRITFI6.231 
00  20  K>1  ,NUM 
HRITEI6.21 IN.XINI.YINI 

FORMAT! 10X.5HNIRE  .I2.2SH  X AND  Y COORDINATES  ARE  ,2F]0.2I 

CONT INUE 

URlTE(bt23l 

FORMAT!*  *J 

00  9 M*1  , NUM 

DO  <•  K 1 1 *NUM 

IF IK-Mllb.17,16 

Z IK,M|=bO.»ALOGI I I XIK1-X IF i>»«2« ( Y|K l ♦YIMI I *«21»»0. 5/ II X|K|-XIM1 i 
l*«2* IYIK!-Y1M)I*«2I»«0.5) 


00  TO  9 

ZIK,M|*bO.«ALOGI2»V{Kl/RADiKl  1 

CONT  INUE 

DO  5 M -1 ,NUM 

00  5 K * 1 , NUM 

L*L*  1 

COIFl(L)*ZIK,MI 
COM  INUE 

DO  b K 1 1 tNUM 
C0EFR1K ) = Y|« I 
CONT INUE 
EPS*] .0E-07 

CALL  01  LG  I COE FR .COEFL.NUM, l, EPS, IER I 
DO  7 N*  1 » NUM 
CUR  I N I *COE  FR IN  I 
CONT INUE 

-4)0^28  M*  1^ N JM  •»  — • -* 


DO  27  N * I ,NJM 

VOLT IM  I *Z  1 M »N  1 »CUR  IN  I « VO  IT  |MI 
CONT  INUE 
CONTINUE 
DO  29  N *1  * NJM 

ABCIN)*ZIN,N l«CUR IN) /VOLT  INI 

HRITEI6.30IN.ABCINI 

CONT  INUE 

F0RMATI10X.36HLEVEL  OF  ARRAY  INTERACTION  FOR  WIRE  .12, 3H  * .E10.91 
HR  IT  E lb  »2  3 ) 

CURR  *CUR ill 

DO  2 5 N * 1 , NUM 

IF ICURR-CURINll2b.2S.25 

CURR  *CUR I N I 

CONT  INUE 

DC  9 N = 1 » NUM 

ACUR  IN  I *CUR I N I / CURR 

HR  | T E lb >10  IN , AC  UR  INI 

FORMAT  1 10X(20HC UR RENT  ON  CONDUC TOR  . 1 3 , 3H  * .E10.91 

CONT  INUE 

MRITFtb,23l 

DO  11  N*l » NUM 

DO  12  M * 1 t NUM 

ZZ IN I *Z  Z I N I ♦ Z I N ,M  I »C  UR  I M I /C JR  INI 


12 


I - I 

H 


kSSSMHBua 


■ ' .1 ' ■ 


ft 


I 

1?  CONTINUE 

HAITf (6,I31N,Z7(NI 

13  F CR  H A T I 1 0 X ,3?H  IMPEDANCE  TC  GROUND  CF  C ON  0 UC  T OR  . I 3 , 3H  * ,E10.N) 

11  CONTINUE 

*RITEI6,?3) 

DO  1<.  N - 1 ,Nj*< 

ZARR=7ARA«I  ?!(N1 

IN  continue 

ZF  *1  ./Z  AR  R 
«RITE(6,]SIZE 

lb  F0RMAT(10X,30H1MPEDANCE  OF  ARRAY  TC  GROUND*  .ElO.NI 
STOP 
END 


5.  CONCLUSION 


The  wire  array  over  a ground  plan  can  be  handled  in  a fashion 
analogous  to  the  treatment  given  the  shielded  multiwire  cable.  A 
comparison  of  theory  with  experiment  for  a lossy  ground  plane  shows  a 
level  of  correlation  comparable  to  that  observed  with  the  multiwire 
cables  and  adequate  to  justify  the  approach  taken.  It  becomes  a 
straightforward  task  to  extend  the  present  development  to  a code 
adequate  to  meet  the  full  needs  of  an  EMP  vulnerability  assessment  and 
hardening  program. 
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